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We present theories of ‘Natural Neutrinos’ in which neutral fermionic top partner fields are simul¬ 
taneously the right-handed neutrinos (RHN), linking seemingly disparate aspects of the Standard 
Model structure: a) The RHN top partners are responsible for the observed small neutrino masses, 

b) They help ameliorate the tuning in the weak scale and address the little hierarchy problem, and 

c) The factor of 3 arising from Nc in the top-loop Higgs mass corrections is countered by a factor 
3 from the number of vector-like generations of RHN. The RHN top partners may arise in pseudo- 
Nambu-Goldstone-Boson (pNGB) Higgs models such as the Twin Higgs, as well as more general 
Composite, Little, and Orbifold Higgs scenarios, and three simple example models are presented. 
This framework firmly predicts a TeV-scale seesaw, as the RHN masses are bounded to be below 
the TeV scale by naturalness. The generation of light neutrino masses relies on a collective breaking 
of lepton number, allowing for comparatively large neutrino Yukawa couplings and a rich associated 
phenomenology. The structure of the neutrino mass mechanism realizes in certain limits the Inverse 
or Linear classes of seesaw. Natural Neutrino models are testable at a variety of current and future 
experiments, particularly in tests of lepton universality, searches for lepton flavor violation, and 
precision electroweak and Higgs coupling measurements possible at high energy e''"e“ and hadron 
colliders. 


I. INTRODUCTION 

The first run of the Large Hadron Collider has placed 
significant pressure on the hypothesis of natural elec¬ 
troweak symmetry breaking. This pressure comes from 
two sources. First, the measured properties of the Higgs 
boson are consistent with the predictions of the Standard 
Model (SM), whereas one would generically expect devia¬ 
tions from these predictions for a natural Higgs. Second, 
direct searches for new states responsible for softening 
the sensitivity of the Higgs mass to high scales, such as, 
e.g., top partners, have so far turned up empty. In par¬ 
ticular, for the familiar case of colored top partners, such 
as scalar top squarks in supersymmetry or fermionic top 
partners in composite Higgs theories, strong constraints 
have resulted from these searches, in some cases pushing 
their allowed masses into unnatural territory [IH3]- 

In light of this situation, there has been renewed inter¬ 
est in the idea of neutral naturalness, which hypothesizes 
that the partner states responsible for the cancellation of 
the quadratic divergences do not carry color charge. Such 
color-neutral partners are more difficult to constrain di¬ 
rectly at the LHC due to their significantly smaller pro¬ 
duction cross sections, thus allowing more natural the¬ 
ories of electroweak symmetry breaking at the price of 
additional model complexity. This novel approach to nat¬ 
uralness dates back to the Twin Higgs model [4], and a 
number of interesting extensions and generalizations have 
been proposed in the literature [MSI]- Recent studies 
of the experimental constraints and signatures are pre¬ 
sented in Refs. [I71II9]. Color-neutral top partners are 
also increasingly motivated by the observed properties of 
the Higgs boson, which indicate that the Higgs-gluon- 
gluon and Higgs-photon-photon couplings are approxi¬ 
mately SM-like, suggesting the Higgs may not be coupled 
to light colored or charged fields. 


Intriguingly, there are other strong hints in nature for 
new neutral states beyond those present in the SM. In¬ 
deed, two of the most compelling empirical suggestions 
of new physics come from the need to generate neutrino 
masses and from the disparate gravitational phenomena 
pointing towards dark matter. Neutrino masses can be 
elegantly explained by the introduction of new neutral 
fermions - the right-handed neutrinos (RHN) - which 
mix with the left-handed neutrinos via Yukawa interac¬ 
tions P^26j . Furthermore, the simplest dark matter 
candidates consist of new cosmologically stable neutral 
particles. A natural question to ask is whether these 
new neutral states, required to understand these empiri¬ 
cal mysteries, can also play a role in naturalness. Stated 
more simply, can these neutral states be top partners? 
Remarkably, the potential connection between neutral 
top partners and dark matter has already been explored 
in an early paper by Poland and Thaler m , who showed 
that neutral top partners could indeed serve as viable 
dark matter candidates. 

In this paper we focus on the previously unex¬ 
plored possibility that the RHNs are simultaneously 
the top partners responsible for canceling the dominant 
quadratic divergence to the Higgs mass. Since the RHNs 
are fermions, we are naturally led to consider theories 
in which the Higgs is a pseudo-Nambu-Goldstone bo¬ 
son (pNGB) [271 - IH^ . in analogy with the pions of QCD. 
Such theories allow for fermionic top partners, which are 
united with the top in a multiplet transforming under 
the spontaneously broken global symmetry of which the 
Higgs is a low energy remnant. In order for the top part¬ 
ners to be neutral under color, the SU(3)c factor of the 
global symmetry must be enlarged to contain, in the most 
straightforward cases, an additional SU(3) factor, which 
is needed to account for the multiplicity factor in the 
top-partner loop. Here we wish to speculate that this is 
in fact the flavor symmetry SU(3)jv which acts on the 
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RHNs, thus making a tentative connection between the 
requirement of three top-partner fields for naturalness 
and the existence of three generations of RHN. Addition¬ 
ally, to enforce the cancellation of quadratic divergences, 
a Zi interchange symmetry or larger SU(6) symmetry 
(which contains SU(3)c x SU(3)Ar) is required. 

The Yukawa couplings of the RHNs to the SM lep¬ 
ton doublets, which are required to generate neutrino 
masses, explicitly break the SU(3)Ar symmetry. There¬ 
fore, unlike other constructions such as the Twin Higgs, 
SU(3)Ar clearly cannot be an unbroken gauge symmetry 
at low energies. Instead, to be consistent with neutrino 
masses, SU(3)Ar can be either a spontaneously broken 
gauge symmetry or an explicitly broken global symme¬ 
try. Since SU(3)Ar is broken while SU(3)c is unbroken 
(and gauged), quadratic divergences will appear at two 
loops, potentially threatening the naturalness of this sce¬ 
nario. However, as we will argue below, these quadratic 
divergences lead to tunings which are tolerable, 0(10%), 
and can even be reduced by gauging the SU(3 )at symme¬ 
try and spontaneously breaking this symmetry at a scale 
close to the RHN mass. 

Regarding the neutrino sector, once we assume that 
SU(3)Ar is broken, various Yukawa interactions and Ma- 
jorana mass terms can be present in the Lagrangian. The 
Majorana mass terms softly break the global symmetry 
and do not introduce new quadratic divergences. The 
neutrino Yukawa couplings on the other hand represent 
a hard breaking of the global symmetry and will lead to 
quadratic divergences. However, provided these Yukawa 
couplings are not too large, these contributions will also 
be under control. This is analogous to the other light 
fermion Yukawa couplings, such as the bottom quarks, 
which do not jeopardize the naturalness of the theory. 
Motivated by these considerations, and the goal of gener¬ 
ality, in this work we will allow all renormalizable terms, 
such as masses or Yukawa couplings, which explicitly 
break the required global symmetry but in such a way 
that the breaking does not spoil the naturalness of the 
theory any more than is already present due to the light 
fermion Yukawas. The approach taken here is thus a phe¬ 
nomenological bottom up exploration of the general class 
of models. 


The general features and bottom up requirements on 
the Natural Neutrinos scenario relating to naturalness 
are outlined in Sec. m In Sec. m we discuss the cou¬ 
plings necessary for the generation of neutrino masses, 
and in Sec. |IV| we describe several explicit low energy 
coset models which motivate the overall structure of the 
models. Although we aim to focus on the phenomenol¬ 
ogy and low energy structure of this framework, we also 
comment on the embedding of these bottom up models 
into standard frameworks including composite and Twin 
Higgs scenarios. In Sec. |V] we will survey various as¬ 
pects of the phenomenology of Natural Neutrinos beyond 
the generation of neutrino masses, including current con¬ 
straints and future experimental prospects. We conclude 
in Sec. VI An Appendix provides a detailed discussion 



FIG. 1: One-loop quadratically divergent corrections to the 
Higgs mass from top quarks and RHNs. These contributions 
cancel in the proposed models as the coupling relation is en¬ 
forced by a global symmetry. 


of proton decay and baryon number conservation in our 
scenarios. 


II. RIGHT HANDED NEUTRINOS AS TOP 
PARTNERS 

We begin by describing the basic ingredients needed 
for the RHN to cancel the dominant top quark-induced 
quadratic divergence^ of the Higgs mass, i.e., to be a top 
partner. This is demonstrated schematically in Fig. 
Given that the RHNs are fermions, we are led to con¬ 
sider theories in which the Higgs arises as a pNGB, which 
provides a framework in which fermions can serve as top 
partners. In such theories, a large global symmetry G is 
broken to a subgroup H at a scale / ^ TeV, and the 
Higgs boson is identified as a pNGB of the coset G/H. 
The essential observation is that the RHN and the top 
quark can be joined in a multiplet transforming under G, 
such that, up to additional symmetry breaking effects, 
the top Yukawa coupling respects the global symmetry. 
In this way, the radiatively induced Higgs mass is gov¬ 
erned by the mass of the RHN top partner, tum ~ A*/ < 
TeV, rather than the ultraviolet (UV) cutoff A ~ 5 — 10 
TeV. 

Let us expand on this idea. The low energy dynamics 
of the pNGBs are described by a non-linear sigma model, 
with the field 


S = e™/-''Eo, n = 7r“T“, (1) 

where 7r“ are the pNGBs, T“ are the broken genera¬ 
tors of G, and (E) = |Eo| = /. The minimal require¬ 
ments on the symmetry breaking pattern are that the 
unbroken subgroup H contains the electroweak group, 
SU(2 )wxU(1)f, and that there are at least four pNGBs, 
dimG — dimiJ > 4 in order to furnish a scalar elec¬ 
troweak doublet at low enegies. We will describe explicit 
coset constructions in Section IlYl 

The third generation weak doublet quark q and weak 
singlet top quark will be embedded in multiplets of 


^ In a full UV completion the Higgs mass is calculable, and the sen¬ 
sitivity to the cutoff A is replaced by a sensitivity to the physical 
threshold, i.e., the masses of the states in the UV theory. 
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G which we denote Q and Q‘^, respectively. In addi¬ 
tion to the top quarks, these multiplets contain neu¬ 
tral top partners N and which we will identify 
with the RHNs. For this to occur, the SU(3)c factor 
of the global symmetry must be enlarged in order to 
accommodate neutral states. The simplest choices are 
SU( 6 ) D SU(3)c X SU(3)jv or SU(3)c x SU(3)Ar with a Z 2 
interchange symmetry, where we have identified SU(3)Ar 
as the flavor symmetry of the right handed neutrinos. It 
is instructive to write a “simplified model” for the top 
quark Yukawa coupling, which can be realized in explicit 
G/H cosets (see Sec. IV): 

( 2 ) 


C = AtQSQ°-kh.c. 


— A+ 


q-ht\ + / ( 1 - p ) N^Nt + . . . 


-I- h.c.. 


where in the second line we have expanded out the G 
multiplets Q, and E in their component fields to 
0{h^h). In Eq. ([^, the SU(3)c index A = 1,2,3 and 


the SU(3)Ar index i = 1,2,3. As we will see in Sec. IV 


in concrete realizations the precise structure of the inter¬ 
actions above can be generalized to include interactions 
with SU(2)w X U(l)y charged states, and with different 
values of the coefficients in front of the couplings, but 
this simplified example will serve to illustrate the basic 
features relating to naturalness in this framework. 

With the additional coupling of the Higgs to the RHNs 
in Eq. ([^, the radiative contribution to the Higgs mass 
parameter, /i^, is softened from a quadratic to logarith¬ 
mic sensitivity. 




87r2 
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This is also depicted in Fig. A naive estimate of the 
tuning in this theory is given by \2 6^‘^, which is 
of order 10% for / ~ 700 GeV and A - 5 TeV. Without 
the RHN top partners, the radiative correction to the 
Higgs mass parameter is quadratically sensitive to the 
UV cutoff, leading to a tuning at the sub-percent level 
for the same choices of / and A. In Eq. (§ and below 
we evaluate the dimensionless couplings at the UV scale 
A as is suggested by a renormalization group-improved 
analysis; see the discussion in Ref. m- 

To generate neutrino masses we must ultimately break 
the SU(3)Ar symmetry. The implementation and break¬ 
ing of the SU(3)Ar symmetry is also relevant for tun¬ 
ing, as has been emphasized recently in Ref. [H]. If we 
take the SU( 3 ) 7 v symmetry as a global symmetry then 
the quadratically divergent two-loop corrections to the 
Higgs mass coming from loops of top quarks which in¬ 
clude SU(3)c gluon exchange are not cancelled. This cor¬ 
rection is given by 


6fj,^ 


3Atg3 .2 
87r4 


(4) 


which, for a UV scale A ~ 5 TeV, yields a tuning of order 
10%. Beyond this correction, parametrically similar two- 
loop contributions arise from the fact that the Yukawa 


couplings to top quarks and top partners run differently 
due to the absence of SU(3)Ar gluons. In [19] it was esti¬ 
mated that the case of a global SU(3)jv symmetry would 
correspond to a tuning of 0(10%) for A ~ 5 TeV. Thus 
in our basic bottom-up approach for now we will sim¬ 
ply assume that SU(3)Ar is a global symmetry and allow 
this level of tuning. We will also allow explicit SU(3)jv 
symmetry breaking terms which lead to the generation 
of neutrino masses, restricting the sizes of such terms in 
such a way that does not increase the tuning further. 

We consider a tuning of 0(10%) to be acceptable for 
the purposes of this work and in the following sections 
focus on the global symmetry limit. However, we note 
that it is possible to reduce this tuning further by instead 
considering the case of a spontaneously broken SU(3)Ar 
gauge symmetry, with gauge coupling (A) « 53 (A). 
In this case, the two-loop top-gluon corrections will be 
largely compensated by symmetry-related contributions 
coming from the top partners and SU(3)Ar gluons, soft¬ 
ening the UV sensitivity of the Higgs mass parameter to 
the scale of SU(3)Ar breaking. If this scale is not too 
far from the G ^ H scale / ^ TeV, the tuning will 
be appreciably reduced. This approach brings with it 
additional model-building questions related to the sector 
responsible for the breaking of the SU(3)jv gauge symme¬ 
try, which, while interesting, will not be pursued further 
in this work. 


There are a number of additional considerations re¬ 
quired for a viable, natural theory of electroweak symme¬ 
try breaking. Besides the correction from the top loop, 
there is an additional important contribution to the Higgs 
mass parameter from the weak gauge boson loops, which 
can dominate over that of Eq. ^ depending on the scales 
/ and A. This loop can be regulated with additional 
gauge boson partners, as occurs in Little Higgs or Twin 
Higgs theories. However, whether or not such gauge part¬ 
ners are required depends again to some degree on the 
amount of tuning one is willing to accept, and in Section 


IV we will consider example models both with and with¬ 


out gauge boson partners. There is also the question of 
generating an appropriate scalar potential for the Higgs 
doublet, h. This occurs due to the explicit breaking of 
the global symmetry by the gauge and Yukawa interac¬ 
tions, but depends in detail on the low energy spectrum 
and to some extent the UV completion. As these issues 
are model dependent, we will set them aside for now and 
move next to the generation of neutrino masses in our 
scenario. 


III. NEUTRINO MASSES 

Through the top Yukawa interaction in Eq. (§ the 
RHN top partner fields N, obtain SU(3)Ar-symmetric 
vector-like masses ~ ytf- However once SU(3)Ar and 
SU(3 )l (the flavor symmetry associated with the three 
SU( 2 )w lepton doublets) are broken, which they must 
be in order to generate neutrino masses and mixings. 
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a number of additional terms may arise consistent with 
the broken symmetries, including the spontaneously bro¬ 
ken electroweak symmetry. There are possible Majorana 
mass matrices from breaking of SU(3)Ar, 

^ ^ ^ , (5) 

and there are Dirac masses which require the breaking of 
SU(3)Ar, SU(3 )l, and electroweak symmetry. These are 

C D + MD,^JN,VJ . (6) 

Combining all terms, we may write the neutrino mass 
matrix in the {v, N, N‘^) basis as 




Md 

Mh \ 

M = 

Ml 

Mm 

Mm 


\m^J 

Ml 

M^mJ 


where each entry is a 3 x 3 matrix. Since the vector-like 
RHN mass Mjv is SU(3)Ar symmetric, while the other 
terms arise only after SU(3 )n symmetry breaking, we ex¬ 
pect those terms to be suppressed in comparison to Mjv. 
Furthermore, the mass Mn ^ Xtf is predicted to be be¬ 
low the TeV-scale by naturalness arguments. Therefore, 
the hypothesis of RHN top partners robustly predicts a 
TeV-scale seesaw mechanism! 

We now explore the conditions required for the gener¬ 
ation of the light neutrino masses. The determinant of 
the mass matrix Q is given by [3S] 

\M\ = \Mn\x\MJ;-M2jM^^Mm\ (8) 

+ {Md - M'oM^^Mm) 

where standard matrix multiplication is understood 
within the individual determinants. This can be under¬ 
stood schematically from the determinant in the case of 
a single generation, 

\M\ ~ MoiMl.MD-MNM't) (9) 

- MnMd) . 


Alternatively, one can integrate out the heavy RHNs, 
which leads to the following mass matrix for the light 
neutrinos 


M, 


{Md Mf,) 


Mm Mm \ ^ ( M^ \ 
M], Ml,) 1 , ) ■ 


( 10 ) 

It is clear that it is a collective symmetry breaking which 
leads to the generation of non-zero light neutrino masses. 
In order to generate neutrino masses at least one of the 
Dirac mass terms of Eq. § must be non-zero and it 
is necessary that terms involving both and N are 
present in order that lepton number is broken collectively. 
The minimal options for generating neutrino masses are 
to have non-vanishing entries in the pairs of matrices 


shown in Table [Tj and we also display the approximate 
value of the light neutrino masses generated. Alternative 
options such as non-vanishing {Mm,M^), {Md,Mm), or 
{Mfj,M^) are not sufficient. The former case is obvious 
as there is no coupling to the left-handed neutrinos. The 
latter cases arise essentially due to a rank condition. It is 
interesting to note that the first option in Table |l] is sim¬ 
ilar to the so-called ‘Linear’ [30] seesaw while the second 
and third option provide examples of the ‘Inverse ’ m 
seesaw (see also Refs. [38ll4T] for variant Inverse see saw 
scenarios). One particularly novel aspect of the collec¬ 
tive breaking of lepton number is that it allows for com¬ 
paratively large neutrino Yukawa couplings, resulting in 
various exotic phenomenological consequences. We will 
explore this feature in detail in Section [Vj 


TABLE I: Minimum pairs of mass matrices required for neu¬ 
trino mass generation. The Linear seesaw [36] may be realized 
for a pair of Dirac mass matrices and the Inverse seesaw m 
for one Dirac and one Majorana mass matrix. 


Non-zero mass-terms 

Approximate neutrino masses 

MD,M't 



~ MlMulMl 

Mm, Mb 



Although this work is not concerned with the details 
of flavor structures in the SM, this class of neutrino mass 
models may also be interesting from the perspective of 
masses and mixing angles. Large hierarchies of masses 
are observed in the charged lepton and quark sector and 
also in the quark sector there are apparent hierarchical 
structures in the mixing angles. However, in the mod¬ 
els presented in this work the neutrino mass eigenvalues 
and mixing angles arise as the product of at least two 
seemingly unrelated matrices, suggesting that the mass 
and flavor structure of the neutrino sector would be very 
different from the other fermions of the standard model. 

Given that the neutrino mass terms in Eqs. if,® 
break the SU(3) n symmetry, it is reasonable to imag¬ 
ine that they originate from non-renormalizable opera¬ 
tors. In this case, an immediate question is why we do 
not simply write down the Weinberg operator (L • h)^/A 
[42] . Assuming that the couplings in Eqs. are the 
only spurions breaking SU(3)Ar and U(1 )l, then clearly 
any light neutrino mass must be proportional to the spe¬ 
cific products of the spurions listed in Table (|T]). However, 
the scale suppressing these spurions may be much higher 
than niM, cor responding to the UV dynamics that gener¬ 
ate Eqs. ®, justifying our neglect of a bare Weinberg 
operator. 

In specific models, there may be further symmetries 
which constrain the form of the neutrino mass terms in 
Eqs. 0,@. Notably, in models with gauge boson part¬ 
ners a la Twin Higgs, there are additional gauge symme¬ 
tries under which the top partners are charged. Thus, to 
generate some of the terms in Eqs. requires addi¬ 

tional insertions of the fields which spontaneously break 
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this symmetry. From the bottom-up approach we are 
pursuing here, there is no obstacle in writing down such 
terms. 

The Dirac masses Mo , originate from Yukawa cou¬ 
plings of the form y^LkN, etc. and thus represent a hard 
breaking of the global symmetry protecting the Higgs 
mass. This is also true of the other SM fermions, pro¬ 
vided we do not embed them in a G multiplet with their 
own partner fields. However, this hard breaking does not 
upset the naturalness of the theory provided the Yukawa 
couplings are small enough. The radiative correction to 
the Higgs mass parameter from these Yukawa couplings 
is given by 




_Ma 2 

87r2 ■ 


( 11 ) 


For a UV cutoff A ^ 5 TeV, the tuning estimate is greater 
than 0(10%) for neutrino Yukawa couplings smaller than 
ijv < 0.25. The Majorana masses Mm,M^, softly break 
the global symmetry G, and thus do not introduce new 
quadratic divergences. However, the Higgs mass parame¬ 
ter is still sensitive to these parameters through the phys¬ 
ical mass of the RHN top partners, and therefore natu¬ 
ralness requires Mm,MI^ < TeV. We will impose these 
constraints on the size of the explicit SU(3)Ar breaking 
parameters generating neutrino masses throughout this 
work. 


IV. MODELS 


a number of occasions [331 HIST]. As a scenario for neu¬ 
tral top partners, this model was proposed in Ref. m, 
and we will largely follow their discussion here. The 
global symmetry is taken to be G = SU(6) x SU(3)w- 
The SU(6) factor contains SU(3)c, which is gauged, and 
the SU(3)Ar flavor symmetry of the RHN. The SU(3)w 
factor is broken to the SU(2)w subgroup, resulting in 5 
pNGBs, 4 of which form an SU(2)u/ doublet to be identi¬ 
fied with the Higgs. We also note there is a gauge singlet 
pNGB, although it will not be important for our discus¬ 
sion. 

The global symmetry breaking is induced by a scalar 
field, E, transforming as a 3 under SU(3)w, which ac¬ 
quires a vacuum expectation value Eq = (0,0,/). The 
pNGBs can be parameterized by the non-linear sigma 
field as in Eq. Q, with 

/ 0 0 hi \ 

n = 7r“T“ = 0 0 ha +..., (12) 

\h\ hi 0 J 


with T“ the broken generators of SU(3)vi/ and we have 
suppressed the singlet pNGB. We may write the sigma 
field explicitly as 


E = 


/ 

V 


ih- 


ih' 


sH\h\/f) 

\h\/f 

sin(|h|//) 

\h\/f 

/cos(|h|//) 


\ 

/ 


(13) 


With the basic framework for a RHN top partner set 
out, we now turn to explicit phenomenological models 
where the low energy dynamics of the Higgs and the 
top partners are embedded within a specific symmetry¬ 
breaking structure. We will focus on three minimal sce¬ 
narios for a pNGB Higgs: a) SU(3)/SU(2), b) the cus- 
todially symmetric SO(5)/SO(4), and c) SU(4)/SU(3). 
Our main purpose in this section is to exhibit possible 
embeddings of the RHNs within the top quark multi- 
plets, which will enforce an interaction structure analo¬ 
gous to Eq. ([^ leading to the cancellation of quadratic 
divergences. We will not wed ourselves to any particular 
UV physics, with the hope that these models may find 
both strongly coupled completions {e.g. composite Higgs) 
or perturbative completions {e.g. SUSY). Furthermore, a 
detailed investigation of radiative electroweak symmetry 
breaking is beyond the scope of this work. We delay a 
discussion of the phenomenological constraints related to 
the EWSB sector, such as Higgs couplings measurements 
and electroweak precision tests, as well as observable sig¬ 
natures to Sec. El 


A. SU(3)/SU(2) 

The minimal coset that furnishes a doublet pNGB un¬ 
der SU(2)vy is SU(3)/SU(2) and has been investigated on 


where \h\ = VhJh. 

We now introduce the top quark and their partners, 
the RHNs. Following Ref. [II] we add two fields Q, Q‘^, 
which transform under SU(6) x SU(3)iy as Q ^ (6,3), 
Q'^ ~ (6,1). We write these multiplets as 



The top Yukawa coupling is written in an SU(6) x SU(3) w 
symmetric manner as 


C = AtQEQ" 


= At 


iq^hfX + / 



WiV/ -h ... 


(15) 
-1- h.c.. 


which precisely reproduces the structure of our “simpli¬ 
fied model” in Eq. (Ill- 

Notice that Q is an incomplete multiplet under the 
global symmetry group G. Naively this is worrisome since 
for generic incomplete multiplets the global symmetry is 
explicitly broken, and the cancellation of quadratic di¬ 
vergences does not hold. However, in this case one can 
understand the cancellation as a result of “twisting” m- 
starting from a full multiplet under SU(3)c x SU(3)w, we 
have twisted one component, V, from SU(3)c to SU(3)Ar. 
In this way, the couplings are structured such that the 
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pNGB Higgs mass is protected against quadratic diver¬ 
gences. One possible realization of such incomplete mul- 
tiplets is from a compactified 5th dimension, in which 
zero entries in the multiplets in Eq. (14) corre¬ 

spond to would-be zero modes projected out by boundary 
conditions. 


Hypercharge can be accommodated in a straight¬ 
forward manner by enlarging the global symmetry to 
U(6) X U(3), such that Y = with Yq = 

diag(|, |,0,0, 0) is a U(6) generator not contained in 
the SU(3)c subgroup and Ya = diag(|, 0) is a U(3) gen¬ 
erator not contained in the SU( 2 )i 4 / subgroup. Finally, 
it is clear that this model, while being minimal, does not 
afford the possibility of having gauge boson partners and 
therefore the quadratic divergence coming from SU(2)w 
gauge bosons is not cancelled. 


B. SO(5)/SO(4) 


In this model the symmetry breaking pattern is 
SO(5)/SO(4), which furnishes exactly four pNGBs which 
transform as a (2,2) under the unbroken SU( 2 )i x 
SU(2)i^ ~ SO(4) symmetry and are identified with the 
Higgs [31]. In the context of composite Higgs theories, 
this scenario is attractive due to the custodial symmetry 
of the strong sector, which protects against large contri¬ 
butions to the p parameter. As before, we also assume an 
SU(6) global symmetry containing SU(3)c and SU(3)Ar. 
Note that in this scenario, as in the previous one, there 
is no possibility of neutral gauge boson partners, in con¬ 
trast to the recent proposals of Refs. nannum based 
on the coset SO(8)/SO(7). See also Ref. jlH] for pNGB 
composite Higgs models exploring the interplay between 
top partners, the lepton sector, and naturalness. 

The global symmetry is broken by a scalar field, E, in 
the 5 representation of S'0(5), which acquires a vacuum 
expectation value Eq = (0, 0, 0, 0, /). The Goldstones are 
parameterized as 

E = Eoe"®/^, H = y2/i“r“, (16) 


where T® are the broken generators. Explicit expressions 
for the generators are found in Ref. jUj. Using these 
generators, one obtains the expression 


^ sin|fe|// 
\h\lf 


/ 

V 


hi 

h2 

hs 

hi 

\h\cot\h\/f 




/ 


(17) 


where \h\ = = y/2h^h 

We now consider the top sector. We embed the left 
handed top and bottom into Q ~ ( 6 , 5) and the right- 
handed top into ~ (6,1). Each of these multiplets 


also contain top partners. The explicit embedding is 


Q 


1 

71 


b —ib t 
£ i£ N 


it 0 
-iAf V2N 






(18) 

where a complete SO(5) multiplet has been “twisted” be¬ 
tween the upper and lower rows in Q, in analogy with the 
SU(3)/SU(2) model described above. Therefore, while Q 
is an incomplete multiplet under the full global symme¬ 
try, the top quark Yukawa and top partner-Higgs cou¬ 
plings are structured so that the pNGB Higgs mass is 
protected. 

The top Yukawa is written as 


C D AtQEQ^-hh.c. 


(19) 


D At [q^hY^ + N^N^] + hx. 

In comparison to the model of Sec. |IV A| based on 
SU(3)/SU(2), this model contains an additional elec- 
troweak doublet of fermions, L ~ (1,2, i) = {£,N) (to 
be distinguished from the SM lepton doublet L). 

It is interesting to examine how the quadratic diver¬ 
gences cancel in this model. The second term in Eq. (20) 


yields the same contribution to the quadratic divergence 
as the top quark loop. However, the coupling of the 
term is twice as large as the one in Eqs. ([^ and 
( |24| ), and so the NN'^ loop neatly cancels the quadratic 
divergences. 

As it stands, some of the fermions are massless, but 
we can add mass terms to the Lagrangian, such as, e.g., 
where = (A/'“,£l°). These mass terms softly 
break the global symmetry, but naturalness will not be 
spoiled provided that these masses are not significantly 
larger than /. Note also that there are additional pos¬ 
sibilities for neutrino mass terms in this model in com¬ 
parison to the “simplified model”, although the discus¬ 
sion there can be straightforwardly generalized to include 
these new terms. 

In this scenario, hypercharge can be realized as follows. 
The global symmetry is enlarged to U( 6 ) x SO(5), such 
that Y = Y^+Tf^, with Yg = diag(|, |, |, 0, 0, 0) is a U( 6 ) 
generator not contained in the SU(3)c subgroup and 
is the generator associated with the SU(2)fl subgroup of 
SO(5) (see Ref. [33] for the explicit form of this genera¬ 
tor). As in the SU(3)/SU(2) model, this model does not 
contain gauge boson partners. 


C. SU(4)/SU(3) 

Finally we discuss a model based on the Twin Higgs |4], 
which was the first model to realize neutral top part¬ 
ners. The global symmetry is taken to be SU( 6 ) x SU(4) 
(or alternatively [SU(3) x SU(2)]^ with a Z 2 interchange 
symmetry). The SU(4) symmetry is assumed to be 
broken to SU(3) by the vacuum expectation value of a 
scalar field in the fundamental representation of SU(4), 
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So = (0,0,0,/), yielding 7 pNGBs. The SU(2) x SU(2) 
subgroup of SU(4) is gauged such that 3 of the pNGBs are 
eaten by the Twin SU(2) gauge bosons. The remaining 4 
pNGBs form a doublet under SU(2)vi/ and are identified 
with the Higgs boson. The pNGBs can be parameterized 
as in Eq. Q, with 


n = 7r“T“ 


/ 0 0 0 hA 

0 0 0 /i2 

0 0 0 0 

\h\ hi 0 0 J 


( 20 ) 


with T“ the broken generators of SU(4). We can ob¬ 
tain the expression for S by expanding the exponential 
in Eq. Q 


S = 


^ \h\lf 

sm(|h|//) 

0 

V /cos(|/i|//) J 


( 21 ) 


with \h\ = Vh'^h 

The top quark and its partners are embedded into rep¬ 
resentations of SU(6) X SU(4) as Q ~ (6,4), ~ (6,1). 

These fields are written explicitly as 


mass not only from loops of the light degrees of freedom, 
but also from resonances near the compositeness scale. 
See also the earlier discussion of SO(8)/SO(7) symmetry 
breaking pattern in Refs. [HEllIl]. It would be in¬ 
teresting to consider RHN neutrino top partners in these 
frameworks. 

It should be noted, especially in the Twin Higgs sce¬ 
nario, that there may be other SM-neutral fields which 
could also play the role of RHN. However, here the work¬ 
ing assumption is that only the Twin Top quarks take 
this role. The presence of the additional Twin sector 
fields would also typically lead to significant modifica¬ 
tions of the RHN phenomenology by introducing addi¬ 
tional decay chains and production mechanisms. We will 
comment on these possibilities where appropriate in the 
next section. 


V. PHENOMENOLOGY 

We now consider the main phenomenological conse¬ 
quences of the Natural Neutrinos scenario. Unless specif¬ 
ically noted, throughout we will assume that only the SM 
and RHN top partner fields are present. 


A. Neutrino oscillation data 


Q = 



A/ ) 


( 22 ) 


where L = {E, N) (which are Twin-sector fields to be 
distinguished from the SM lepton doublet L) and N‘^ are 
the neutral top partners. The top Yukawa coupling is 
written in a SU(6) x SU(4) invariant way as 


C = XtQ'EQ^ + h.c. 

= At 


h^h 




(23) 
-I- h.c., 


which provides another example that reproduces the 
“simplified model” in Eq. ([^, ensuring the cancellation 
of quadratic divergences in the top sector of the theory. 

In the usual implementation of the Twin Higgs, the 
second SU(3) Twin color symmetry (here identified with 
SU(3)Ar) is gauged, and left unbroken. However, as we 
have described in detail in Sec. |Hj we take this to be 
a global symmetry. However, we do gauge the twin 
SU(2)w symmetry, and therefore the model predicts 
gauge boson partners which cancel the quadratic diver¬ 
gences coming from the SM weak gauge bosons. 

Generalizations of the original Twin Higgs scenario 
in the context of composite Higgs models or their holo¬ 
graphic duals have recently been constructed [161 |20l |2T] . 
These models are based on the symmetry breaking pat¬ 
tern SO(8)/SO(7), and provide scenarios in which there 
is custodial symmetry and also protection of the Higgs 


We begin by detailing how the neutrino oscillation data 
can be described in our framework. Given that natural¬ 
ness dictates the large scale Mjv ^ Utf ^ TeV while the 
other entries in Q can only result from SU(3)Ar break¬ 
ing, it is natural to expect the latter to be suppressed in 
comparison to Mjv, leading robustly to the prediction of 
a TeV scale seesaw. Eurthermore, one novel aspect of this 
framework is the fact that at least two additional entries 
are needed to generate neutrino masses (see Table |T|, 
indicating that there is a collective breaking of lepton 
number. This allows in principle for some of the Yukawa 
couplings to be relatively large, leading to the possibility 
of additional novel phenomena correlated with neutrino 
mass generation, as we will describe below. 

As already mentioned, there are two interesting limits 
in the mass matrix 0 which realize the Inverse seesaw 
(when or M/,, Mm are the only additional non¬ 

vanishing entries) and the Linear seesaw (when Mjj, M/, 
are non-vanishing and the Majorana masses vanish). For 
concreteness, we will specialize to the case of the Inverse 
seesaw in what follows, with the neutrino mass matrix 

f 0 Md 0 \ 

M = \ 0 Mn \- (24) 

VO < / 

In the limit Mn, the spectrum consists of 

three light neutrinos and three pairs of heavy pseudo- 
Dirac fermions. This can be easily seen in the limit of 










one generation, which upon diagonalization of (24) leads 
to the eigenvalues: 




mN,± 


Ml ’ 

± ( Mjv + ^ 
' 2M. 


N 


M 


M 


(25) 

(26) 


valid in the limit Md, Mpf. We observe the col¬ 

lective breaking manifests through the dependence of rrii, 
on both M]j and The heavy states are split by an 

amount |mAr_+| — |mAr__| « For a natural value 

itin ~ Utf ~ 700 GeV, the light neutrino mass scale of 
order rui, ^0.1 eV can be obtained in the two extremes 
of parameter space: 1) y^, ~ 0.1, ^ 10 keV and 2) 

~ 3 X 10-®, ~ 100 GeV. 

The current global fit to the standard three flavor oscil¬ 
lation scenario for the normal (inverted) ordering yields 
the following ranges for the squared mass differences and 
PMNS mixing angles gaEO]: 

Amil [10-5 eV"] = 7.50;°;[? (7.50l°;[?), 

Ami (Ami^) [IQ-^ eV^] = 2.547l°;°47 (_2.449;n4?). 
sin^ 012 = 0.304lO;O[3 (O.304t0;0[3), 

sin^ 023 = 0.452™ (0.579tnf7), 

sin2 0i3 = O.O218t°:°°[° (0.0219™J), (27) 

while the GP-violating phases are currently uncon¬ 
strained. After integrating out the heavy RHNs we ob¬ 
tain the effective mass matrix for the three light neutri¬ 
nos: 

M, - (28) 


A convenient way to automatically reproduce the low 
energy data is to employ the Casas-Ibarra parameteri¬ 
zation m of the Dirac mass matrix (see also Ref. [52)1. 
which for the case of the Inverse seesaw is 


Md = Mn (M^)-i/2 R (m,)i/2 Ul 


PMNS’ 


(29) 


where R is in general a complex, orthognal matrix satis¬ 
fying R^R = 1, = diag(mj^j, , m,^^) is the physical 

light neutrino mass matrix, and I/pmns is the PMNS ma¬ 
trix. In our numerical results here and below, we will fix 


W;/, DpMNS as allowed by the global data (27), restrict 


ing to the normal ordering for simplicity and scanning 
over the lightest neutrino mass mi,-^. We furthermore 
fix TO AT = 700 GeV while scanning over R and M^. 
In the left panel of Fig. we show the results of the 
scan in the — Mjj plane, where the bar denotes 
the average value of the non-zero entries in the associ¬ 
ated mass matrices. All points in the plot reproduce 
the low energy neutrino data presented in Eq. ( |^ (the 
gray points are excluded by additional constraints to be 
discussed b elow ). The black line represents the correla¬ 
tion of Eq. (25), Md = Mn \Jn^/ notice that the 
points cluster around the line as expected. 


B. Non-unitarity of the PMNS matrix 


A novel aspect of the neutrino mass models proposed 
here, and the Inverse seesaw structure in particular, is 
the possibility of large neutrino Yukawa couplings. Such 
large couplings can manifest as a violation of unitarity 
of the PMNS matrix, which leads to a host of physical 
consequences [55H55] . 

In general, the 9x9 neutrino mass matrix in Eq. 0 
is diagonalized by a unitary matrix, U, such that 

= U'^M U. (30) 

In the presence of the heavy RHNs, the 3x3 sub-matrix 
describing the mixing of the three light neutrinos, U, 
is no longer unitary {U is to be distinguished from the 
standard unitary PMNS matrix C/pmns)- 

Such deviations from unitarity can show up in a num¬ 
ber of measurements. For example, the prediction for the 
Fermi constant, Gp, which is extracted from the precise 
measurement of the muon lifetime, is altered from the 
standard case of a unitary PMNS matrix. The experi¬ 
mentally determined quantity, = 1.1663787(6) x lO”® 
GeV”^ [56], is related to Gf in this case as 

Gl = (31) 

In the standard case of mixing between only 3 light neu¬ 
trinos, PMNS unitarity guarantees {UUlij = Sij, while 
in our scenario this is no longer the case. 

Indeed, a suite of predictions for weak interaction ob¬ 
servables are affected by the deviations from unitarity of 
U, including Z and W boson decays, .Z-pole asymmetries, 
invisible Z-boson width, the W-boson mass, weak mix¬ 
ing angle measurements, lepton flavor universality tests, 
lepton-flavor violating decays, and quark flavor CKM pa¬ 
rameters. A recent study of these effects is presented in 
Ref. |55] . in which bounds are derived on the size of the 
deviations from unitarity on the quantities {UUlij, with 
z,j = 1,2,3. In particular, for the diagonal elements, we 
apply the following conservative 3cr C.L. limits [55] : 

1 - (UUlee < 0.0018, 

1-(GC/^)mm < 0-0007, 

1 - {uulrr < 0.005, (32) 


We note that these limits assume independent variations 
of the elements in the fits, while in our scenario, there 
will be correlations among the non-unitary parameters 
which are expected to change the precise limits. A full 
study of these effects are beyond the scope of this pa¬ 
per, but we expect the bounds in Eq. (32) to be reason¬ 
ably representative of the correlated ones in our scenario. 
Furthermore, there are limits on the off-diagonal (flavor¬ 
changing) elements, which we will return to below when 
we discuss lepton flavor violation. The green points in 
Figure]^ are allowed by the constraints in Eq. (32), while 
the gray points are excluded by a combination of these 
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constraints and those coming from lepton flavor violation, 
as we discuss next. Future possible lepton colliders such 
as the ILC, FCC-ee/TLEP, and CEPC, with their im¬ 
provement in precision electroweak measurements, along 
with a suite of new low-energy experiments testing lepton 
universality, will probe deviations from PMNS unitarity 
at the 10“^ — 10“® level [55] , 

Finally, it would be interesting to investigate further 
the effects of non-unitary mixing on neutrino oscillations, 
see for instance the recent stuy of Ref. [5S] . 


C. Lepton flavor violation 


Another important consequence of the potentially 
large neutrino Yukawa couplings in our neutrino mass 
models are lepton flavor violating processes. A useful and 
comprehensive study of lepton flavor violating decays in 
seesaw models is given in Ref. [SH] (see also Ref. [sniisi]). 
Of particular importance is the ^ > ey branching ra¬ 

tio, which places the strongest constraint on the neutrino 
Yukawa coupling in our scenario. The prediction for this 
branching ratio is [SH] 


Br(/r —>■ ey) 


2567r2 


E 

2=1 


u^^u:,G 


^N,i 




2 


(33) 

where the loop function is G(a:) = —(1/4)(1 —cc) "^[(2x^-1- 
5x^ — x)(l — x) -I- bx^logx]. Here, mN,i denote the 9 
physical neutrino masses, and (b^^) is the muon mass 
(width), which can be found in Ref. [56]. The strongest 
constraint on this branching ratio comes from the MEG 
collaboration. 


Br(/r ey)MEG < 5.7 x 10 (34) 


In the right panel of Figure]^ we display the prediction 
for this branching ratio as a function of the average Dirac 
mass. This measurement places the strongest constraint 
on the Dirac mass, requiring neutrino Yukawa couplings 
Vi' ^ 0 - 1 - 

Looking forward, there are exciting prospects for im¬ 
proving the sensitivity to LFV processes (see, e.o. fS^I^ 
for recent overviews). For ^ —)■ ey, the MEG-II up¬ 
grade will be able to improve the branching ratio limit 
by roughly an order of magnitude [51]. Other prospects 
for LFV involving muons include experiments such as 
COMET, DeeMee, Mu2e, and PRISM, which are ex¬ 
pected to improve the bound on the coherent muon-to- 
electron conversion rate by several orders of magnitude in 
the coming years [521 [53]. Likewise, the Mu3e experiment 
has the capability to improve the bound on the /it —)■ 3e 
branching ratio by four orders of magnitude, probing 
branching ratios down to the level of 10 “^® [521 1^ - 
There are also exciting prospects for r-LFV processes, 
which can be probed by LHCb and flavour factories such 
as Belle-II and Super KEKB [^, as well as LFV Z de¬ 
cays which can be tested at future lepton colliders [64] . 


D. Neutrinoless double beta decay 


A classic signature of the Majorana nature of the neu¬ 
trino is neutrino-less double beta decay. The amplitude 
for this process contains the sum (see e.g. [55H55| L 
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A rr2 ^N,i 


(35) 


where (p^) ~ (100 MeV)^ is the characteristic momentum 
scale of the process. In principle the heavy neutrinos can 
enhance the rate of this process. In particular, neglect¬ 
ing the mixing angles for the moment, we observe from 
Eq. (351 that the contribution of the heavy neutrinos is 


enhanced in comparison to that of the light neutrinos 
by the factor An /~ {p^)/{ ttl^Mn) ~ 10^. Indeed, 
strong bounds have been placed on the size of generic 
mixing elements, |t/eip < I0“®, for TeV-scale RHNs [58] . 

However, since the process violates lepton number, it is 
clear that the rate must be governed by the small collec¬ 
tive breaking of lepton number. This is true even if the 
mixing angles are large, which can happen if is large. 
Working in the flavor basis, it is easy to see that the con¬ 
tribution of the heavy neutrinos requires two Yukawa in¬ 
sertions and one Majorana mass insertion, and therefore 
must be proportional to An ~ ^ 

This can also be seen in the mass basis, in which case 
there is a Hne cancellation between the amplitudes of the 
split pseudo-Dirac states. Therefore, we conclude that 
the contribution from the heavy neutrinos is negligible. 
It is of course still possible that neutrinoless double beta 
decay is observable in this scenario, although it will be 
a consequence of the light neutrinos and in this sense is 
not different from the minimal scenario in which Majo¬ 
rana neutrino masses are described by the dimension 5 
Weinberg operator. 


E. Higgs couplings 


A generic prediction of pNGB Higgs scenarios is the 
modification of the Higgs couplings to SM gauge bosons 
and fermions, which are conventionally parameterized as 


a = 


9hVV 


atvv'' 


c = 


9hff 


(36) 


which are a function the combination . The pre¬ 

cise form of a, c depends on the symmetry breaking pat¬ 
tern and the embedding of the SM fermions in to G 
representations, but are generically of the form a^c ~ 
l — 0{v'^fP). Gurrently, the LHC 7-1-8 TeV data are con¬ 
sistent with the SM values a = c = I and constrain the 
parameter a (c) at the ±10 (40) % level at 3 ct C.L. [TOllTl] . 
Given that the precision on a dominates, we can infer a 
constraint on the symmetry breaking scale / > 400 GeV. 
Looking towards the future, the LHC will eventually be 
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FIG. 2: Neutrino parameter space: Here we display the results of a scan over the Majorana masses and the Dirac masses 
Md via the parameterization of Eq. (291. All points shown reproduce the low energy neutrino oscillation data in Eq. (271, while 
the green points are also allowed by constraints on non-unitarity of the PMNS matrix and lepton flavor violation. In the left 
panel we have displayed the results in the — Md plane (the average Majorana and Dirac masses). The black line indicates 
the correlation expected from Eq. (25l. The right plot shows the prediction for the lepton flavor violating decay /r —>■ ey as a 
function of Md- We also display the current 90% C.L. limit from the MEG collaboration, Br(/i —>■ 67 ) < 5.7 x 10“^® | 57| . 


able to measure a at the few percent level, probing / ^ I 
TeV. Beyond this, future e+e” machines will have the 
capability to measure a at the per-mille level, probing / 
in the multi-TeV range m- 

As the RHN are neutral under the entire SM gauge 
group no additional corrections to the Higgs-gluon-gluon 
or Higgs-photon-photon couplings are expected. Finally 
we note that, in contrast to more general Inverse see saw 
models [73], since the RHN top partners are constrained 
to be heavier than 0(500 GeV) there are no new Higgs 
decays directly into the new neutrino states. 


F. Electroweak precision tests 


A major constraint on the global symmetry breaking 
scale / comes from the precision electroweak data. There 
are several potential contributions to the oblique param¬ 
eters, S and T [74j . First, there is an irreducible contri¬ 
bution which arises due to the modification of the hWW 
and hZZ couplings, which in the models above can be 
described in terms of the parameter a ^ l + 0{v‘^ / de¬ 
fined in Eq. ( |36| ). This modification implies that, in con¬ 
trast to the SM, the IR log divergences in the gauge boson 
vacuum polarizations do not completely cancel up to the 


scale A, leading to a contribution to S and T |7S1|7B|: 




AT = - 




(I-a^)log 


(!)■ 


(37) 


The corrections depend logarithmically on the UV cutoff. 
For instance, in the case of a strongly coupled UV com¬ 
pletion, one expects A ^ 47r/ = 47rri/\/l — a^. Using, 
for example, the most recent results of the Gfitter group, 
S = 0.06 ± 0.09, T = 0.1 ± 0.07, with a correlation coef¬ 
ficient p = 0.91, we find that this places a strong bound 
/ > 1 TeV at 3cr G.L.. However, it should be noted that 
the SM point [S = T = 0) lies outside the Icr ellipse. 
In this light, a more conservative bound can be placed 
by measuring the deviation of AS, AT in Eq. ( jWj ) from 
the central value rather than the SM point, in which case 
the bound of / weakens to / > 500 GeV. If instead one 
considers a perturbative completion the UV cutoff can 
be appreciably lower, reducing the shifts to the oblique 
parameters in Eq. (371. 

Beyond these irreducible pieces, there can be other 
contributions to the oblique parameters depending on 
the UV completion. In strongly coupled UV completions, 
such as composite Higgs theories, one expects tree-level 
contributions from the resonances at the scale rrip ^ Ppf 
where Pp parameterizes the strength of the coupling of the 
resonances |77j. If the unbroken subgroup H does not en¬ 
joy a custodial symmetry, one expects large contributions 
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to the T parameter of order {l/a){v^/p), constraining 
/ to be in the unnatural multi-TeV range. Provided this 
contribution is controlled by a custodial symmetry, one 
still expects a contribution to the S parameter of order 
AS” ~ {Att/ gp{v'^/p), which depending on gp can be 
larger than the contribution in Eq. (37). However, in 


a weakly coupled UV completion, such as SUSY, these 
contributions are expected to be absent. 

Finally, there can be model-dependent contributions 
from other light states in the spectrum, such as the part¬ 
ners to the SM fields. Importantly, one can obtain a pos¬ 
itive contribution to the T parameter if the new states 
break custodial symmetry through their couplings to the 
Higgs, allowing to significantly relax the constraints on 
/• 

Possible future e+e machines, such as the ILC, FCC- 
ee, and CEPC, will be able to determine the oblique pa¬ 
rameters at the percent level [78], probing scales / ^ 1 
TeV. 


Section IV C In particular, in this model the RHN (the 
“Twin top”) may dominantly decay into a Twin W bo¬ 
son and its Twin SU(2) fermion partner E (the “Twin 
bottom”).^ There is a great deal of Twin sector spec¬ 
trum dependence on the possible decay chains. However, 
we note that the Dirac masses of Eq. ([^ may be the 
only couplings that mix Twin sector fermions with SM 
fermions, and thus even if there are long decay chains 
within the Twin sector they may terminate back into the 
visible sector through these operators, potentially lead¬ 
ing to large multiplicity final states involving fermions. 
Furthermore, an additional experimental opportunity is 
present in these scenarios as the Higgs boson can decay to 
these light Twin sector states. In general, the absence of 
a confining SU(3)Ar force suggests a distinct exotic Higgs 
decay phenomenology in comparison to usual Twin Higgs 
scenarios, such as that described recently in Ref. [19] . 


1. Hadron colliders 


G. Direct production at high energy colliders 

It is possible that the RHN top partners may be pro¬ 
duced at colliders. Before considering the possible pro¬ 
duction mechanisms it is worthwhile to first determine 
the decay channels. Here we will consider the minimal 
scenario is which the RHN top partners are the only 
new states in the spectrum. If the Yukawa couplings, 
hNL and hN'^L, were set to zero the left-handed neu¬ 
trino masses would vanish and the RHNs would become 
stable due to a Z 2 symmetry. Thus in a complete model 
all decays must proceed through the Yukawa couplings. 
Making use of Goldstone equivalence the possible decay 
chains of the heavy RHN are 


N, iV“ 
N, iV“ 
N, iV“ 


hu , 

Zv , 

w+i,w-i. 


(38) 


Thus, there are a variety of final states which may be 
useful for collider searches, in particular involving missing 
transverse energy (MET) and/or leptons resulting from 
the fact that decays proceed through lepton Yukawas. 
Furthermore, decays of the heavy RHN will have a final 
state flavor structure which is correlated with the flavor 
structure of the light neutrino masses and mixings. Thus, 
in the charged current decays N ^ Wl final states which 
are not democratic in lepton flavor may arise, giving clues 
as to the origin of neutrino masses and mixings at collider 
experiments. 

We emphasize again that here we are considering a 
minimal spectrum with the only new states being the 
RHN top partners. However, it should be kept in mind 
that if there are additional light neutral states then the 
RHN are likely to decay dominantly into those states. 
This is particularly relevant if the model is embedded 
in a Twin Higgs scenario, such as the one described in 


It is possible to produce the RHN top partners at 
hadron colliders. In the limit of vanishing Yukawa cou¬ 
plings they may only be pair-produced, and this proceeds 
via gluon fusion gg ^ h* ^ due to the coupling 

of Eq. ([^. This production process is relatively model- 
independent as the coupling is dictated by the cancella¬ 
tion of quadratic divergences. Assuming only the RHN 
and the SM the possible final states from RHN pair pro¬ 
duction are 


NN^ 

iViV° 

iViV° 

iViV° 

iViV° 

iViV° 


hhvv, 

hZvv^ 

ZZvv^ 

Wlhv, 

WlZv, 

WlWlv. 


(39) 


In Fig. [^we show the inclusive cross-section for pair pro¬ 
duction at the 14 TeV LHC, which was calculated by 
implementing the model in FeynArts, and then using 
the FeynArts, FormCalc, and LoopTools suite of 
packages izsiiHni to perform the tree-level and loop cal¬ 
culations. These cross sections are typically small, hence 
at most C>(few) RHN pair production events would be ex¬ 
pected with 300 fb~^ of integrated luminosity at the LHC 
for small RHN mass. The pair production cross section 
is significantly larger at 100 TeV suggesting a 100 TeV 
collider would be more promising for discovering RHN 
top partners. 

In the electroweak symmetry breaking vacuum the 
Yukawa couplings mix the RHN with the left-handed 
neutrinos, thus it is possible to singly produce the RHNs 


We would like to thank Gian Giudice for emphasizing this point. 





12 




FIG. 3: The cross section for pair producing the RHN at 14 
and 100 TeV. A vanishing Yukawa coupling has been cho¬ 
sen such that pair production only occurs through the Higgs 
coupling of Eq. §• 


FIG. 4: The cross section for producing a single RHN in asso¬ 
ciation with a charge lepton, pp — ^ Nl, pp — > N‘^1, at 14 and 
100 TeV. A Yukawa coupling of Y = 0.1 has been chosen. 


in association with a v (neutral current) or a charged 
lepton I (charged current). In Fig. we show typical 
cross sections for single RHN production in association 
with a charged lepton, for a relatively large Yukawa cou¬ 
pling tjiy = 0.1.^ For typical RHN masses of 0(700 GeV) 
only a handful of events would be expected with 300 fb”^ 
of integrated luminosity. The cross section increases by 
an order of magnitude when going from 14 to 100 TeV, 
again suggesting a 100 TeV proton-proton would be more 
promising for discovering the RHN top partners. By con¬ 
sidering the possible decay chains of the RHN the typical 
backgrounds to single production in association with a 
neutral or charged lepton would be the di-boson produc¬ 
tion processes 


SIG: Z* -)-vN -)■ hvv 
SIG: Z* uN Zvv 
SIG: Z* ^vN ^ Why 
SIG: W* ^IN ^ hlv 
SIG: W* -)■ IN ^ Zlv 
SIG: W* ^ IN ^ Wll 


BG: hZ{Z vu), 

BG: ZZ{Z ^ lyiy), 

BG: WW(W ^hy), 
BG: hW{W ^hy), 

BG: ZW{W Iv), 
BG: WZ{Z ^ ll). (40) 


By taking into account the typical cross sections for the 
background processes it would appear that detecting the 
RHN top partners may be challenging at the LHG, al¬ 
though it is difficult to assess the possibilities adequately 
without performing a full collider study with optimized 
cuts. A related study in the context of the Inverse see¬ 
saw embedded within the NMSSM finds that it may be 
possible to detect the heavy RHN at the LHC in the 


^ We have compared our calculations against m and find they 
are in good agreement. 


trilepton-|-MET final state with suitable cuts applied [8T] , 
however this was for Mat = 100 GeV which has a much 
larger production cross section. We note that striking 
collider signatures of lepton-number violation are 

suppressed in our scenario by the small collective break¬ 
ing of lepton number. 

Finally, the presence of additional heavy colored states 
in the multi-TeV range may be expected in a UV comple¬ 
tion. Although access to these states at the LHC will be 
challenging it may be possible at a future 100 TeV hadron 
collider, and this possibility deserves further study. 


2. e''"e colliders 


Besides the powerful indirect constraints which could 
be achieved with future e“'"e“ colliders it is worth consid¬ 
ering the direct constraints which may be possible. Due 
to the expectation that u// <C 1 the RHN masses are 
expected to be Mat > 400 GeV, thus pair production of 
RHN would require at least a 1 TeV e+e“ collider such 
as the ILC, and if the masses were Mat > 700 GeV then 
a significant increase in CM energy would be required 
for pair production. On the other hand RHNs could be 
singly produced through e+e“ —)■ Z* —>■ lyN and a 1 TeV 
e“''e“ collider such as the 1 TeV option of the ILC may 
be able to probe a significant portion of the relevant pa¬ 
rameter space. A study for lighter RHN at the ILC was 
performed in [ST] (see also |50|) and it was found that for 
Mat = 150 GeV a statistically significant excess may be 
observable. It would be interesting to perform a similar 
study for the heavier RHN to determine the ty of e.g. the 
ILC to Natural Neutrinos. 
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VI. OUTLOOK 

We now speculate on some aspects that have not been 
considered or developed in this work, but which may lead 
to promising model-building and phenomenological av¬ 
enues of future development. 


A. Neutrino Mass Models 

Our approach to the SU(3)Ar andU(l)L breaking terms 
that generate neutrino masses in this work has been 
bottom-up, allowing all terms at the renormalizable level 
consistent with naturalness constraints. It would be very 
interesting to explore explicit models that generate these 
couplings. Depending on how the SU(3)Ar symmetry is 
implemented, be., as a global or spontaneously broken 
gauge symmetry, this will involve fields responsible for 
the breaking of SU(3)Ar as well as new states mediating 
the couplings between the SM lepton doublets, the pNGB 
Higgs, and the RHN top partners. Explicit constructions 
of this type are likely to give guidance to the expected 
sizes of the effective neutrino Yukawa couplings and Ma- 
jorana masses, and may lead to novel phenomenology. It 
is worth noting that the terms needed for neutrino mass 
generation are small, thus they may be generated by new 
fields well above the UV scale A, or at A but with sup¬ 
pressed couplings. 


B. UV Completion 


The couplings and field content in this work were mo¬ 
tivated by symmetry considerations which ensured the 
cancellation of quadratic divergences due to the IR de¬ 
grees of freedom alone. Ultimately it will be desirable 
to construct a full UV theory with a calculable Higgs 
mass that dynamically accounts for the origin of SU(3)Ar 
and U(1 )l breaking. In such a UV completion, there are 
likely to be additional naturalness considerations that 
go beyond the naive bottom-up coupling structures of 


Sec. 1^ Thus, although the low energy structure of Nat¬ 
ural Neutrino models appears relatively uncomplicated, 
this does not necessarily imply that realizing a fully UV 
complete model would be straightforward and work to¬ 
wards this goal is necessary to put Natural Neutrino mod¬ 
els on a firmer footing. 


C. Leptogenesis 

An attractive possibility would be if the baryon asym¬ 
metry could be explained in the Natural Neutrinos frame¬ 
work from new processes at the weak scale. The models 
may realize the Sahkarov conditions [^. There is ex¬ 
plicit lepton-number violation which feeds into baryon 
number violation while electroweak sphalerons are ac¬ 
tive. It may be possible to achieve sufficient CP-violation 


though the complex phases in the Yukawa couplings 
and/or Majorana mass matrices. The small Yukawa cou¬ 
plings and Majorana masses required for small neutrino 
masses may also lead to out-of-equilibrium processes. In 
fact leptogenesis has been previously found to be possible 
in the Inverse and Linear seesaw models |92l |93] , giving 
support to the possibility that it may be possible in the 
Natural Neutrinos framework. 

There is, however, a potential obstacle. The cou¬ 
pling required by naturalness £ D hh^ NcN/f may keep 
all of the RHN fields in thermal equilibrium in the 
early Universe and this may make it difficult to sat¬ 
isfy one of the Sahkarov conditions. Achieving lepto¬ 
genesis may then require appealing to additional out-of- 
equilibrium dynamics, such as resonant processes involv¬ 
ing the small mass-splittings between the pseudo-Dirac 
RHN. It would be interesting to understand in quanti¬ 
tative detail whether the Natural Neutrino models may 
explain the baryon asymmetry of the Universe. 


D. Dark Matter 

Can some of the top partner fields Ni, Nf be dark 
matter candidates? Naively, it would seem that there 
are no new stable particles since any would-be Z 2 sym¬ 
metry protecting the RHN is violated by the Yukawa 
couplings, allowing the RHN to decay. However, this is 
not necessarily the case. 

For instance, suppose we make N^, odd under Z 2 , 
while all other fields are even. Clearly, the top Yukawa 
coupling in Eq. (§ respects this Z 2 and the Higgs mass 
protection works as before. The Z 2 symmetry will zero 
out entries with only one iVa or field in the neutrino 
mass matrix Q. However, there is in fact still enough 
freedom to generate three light neutrino masses. This 
can be seen by inspecting the determinant in Eq. (|^. 
The expression is written as a product of three determi¬ 
nants of 3 X 3 matrices, and one can see that the first 
two determinants are clearly non-zero due to the pres¬ 
ence of Mjq. The final determinant in Eq. ([^ is also 
non-zero provided the remaining entries in both Md and 
are non-zero. Therefore, in this example, V 3 , V| 
are stable and potential dark matter candidates, while 
simultaneously three light neutrino masses are generated 
through their couplings to the other RHN top partner 
fields. As mentioned in the introduction. Ref. m previ¬ 
ously demonstrated that neutral fermionic top partners 
are potential dark matter candidates, hence this may be 
an interesting future direction of study for the Natural 
Neutrinos framework. 


E. Majorana Top Partners 

One potential area for future development would be 
to realize Majorana RHN as top partners, rather than 
the vector-like RHN arising in the models studied here. 
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In principle it would seem possible to cancel quadratic 
divergences with a term £ D hW /f. However in or¬ 
der to enforce the required coupling at the TeV scale 
it would be necessary to embed the top quarks (which 
are in a complex representation) and the Majorana RHN 
(in a real representation) into an incomplete multiplet of 
some UV symmetry. We did not find such a symmetry 
structure, however it may be possible with further study. 

F. Connecting Nc with Np 

As the Natural Neutrinos scenario enforces equality 
between the number of flavors of RHN and the number 
of colors in QCD {Njq = Nc) it is tempting to specu¬ 
late as to whether it might be possible to tie the number 
of SM fermion families to the number of RHN in some 
encompassing scenario, i.e Np = Njy. Due to the first 
equality such a construction would then realize the very 
attractive possibility that the resolution of the little hier¬ 
archy problem and the mechanism behind the generation 
of neutrino masses leads to the prediction that the num¬ 
ber of fermion families must be equal to the number of 
colors in QCD, Np = Nc- However, we did not find any 
construction that achieves this goal and believe it may 
be difficult, especially as the Natural Neutrinos scenario 
treats the third generation of quarks separately from the 
first two. 


VII. SUMMARY 

In summary, the ‘Natural Neutrinos’ scenario proposed 
here represents a new class of bottom up ‘neutral natural¬ 
ness’ models which address the little hierarchy problem 
while explaining the absence of new colored fields at the 
LHC. In our framework, the neutral top partners are si¬ 
multaneously the RHNs responsible for the generation 
of the light neutrino masses. The models also enforce a 
novel connection between the number of vector-like RHN 
(Abv) and the number of colors in QCD (N^). These mod¬ 
els may arise in a variety of pNGB Higgs scenarios and 
we have sketched three specific models to demonstrate 
this. The precise structure of neutrino mass generation 
is based on a collective breaking of lepton number, al¬ 
lowing for the possibility of large neutrino Yukawa cou¬ 
plings. The models predict a plethora of potential signals 
in low-energy tests of lepton universality and lepton fla¬ 
vor violation, as well as possible signatures at high energy 
collider experiments. 
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Appendix A: On Proton Decay 

In the model descriptions of Sec. |IV| the visible sector 
quarks are described as living within larger multiplets 
(specifically SU( 6 )) that contain also the hidden sector 
quarks. As these hidden sector quarks become the RHN, 
and lepton number is broken, it would seem then that 
the multiplet structure may lead to proton decay as any 
global baryon-number U(1 )b which acts on the full multi¬ 
plets containing visible and hidden sector quarks together 
has clearly been broken.^ In this section we will demon¬ 
strate that proton decay may be avoided in the models 
presented in Sec. |IV[ 

As a warm up let us consider the Twin Higgs model. 
In this scenario the approximate Twin symmetry leads 
to the appearance of a global SU( 6 ) x SU(4) symme¬ 
try, however in reality this is only apparent as there 
are two copies of the SM fields, and the SM quarks 
and Twin quarks do not actually live in SU( 6 ) mul¬ 
tiplets, but rather SU(3) x SU(3). The appearance 
of a full SU( 6 ) symmetry is a consequence of the ex¬ 
change symmetry which identifies the gauge couplings 
of both SU(3) groups as being equal. With proton de¬ 
cay in mind, the (anomalous) global U(l) symmetries 
are U(1 )b x U(1)l x U(I)bt x U(l)iT, where the lat¬ 
ter two denote the Twin symmetries. Thus in order to 
generate the neutrino masses as described the symme¬ 
try is broken to U(1 )b x U(1)l x U(1)bt’ x U(l)iT —>• 
V{1)b X X U(I)iT and Twin baryon number 

breaking does not imply any breaking of the SM U(1 )b 
baryon number and the proton does not decay. 

The situation for the Twin Higgs is similar to the case 
for the other models as in pNGB scenarios typically only 
subgroups of a larger global symmetry group are gauged. 
Thus on a case-by-case basis additional global baryon- 
number symmetries which arise as elements of the full 
symmetry group may typically be imposed which protect 
the proton from decay. However in order to demonstrate 
the origin of proton stability within the context of the 
full symmetry group it will be useful to consider the re¬ 
quired global symmetry structure and discuss the gauged 
subgroups at a later stage. This is useful because once 
subgroups are gauged it is sufficient to demonstrate that 
this gauging does not break the relevant U(l) factors. 
Thus demonstrating a conserved baryon-number in the 
context of the full global symmetry structure is a useful 
first step. 


^ We thank Ben Gripaios for raising this interesting point at the 
CERN-CKC Neutral Naturalness workshop. 
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We will now consider another example with an explicit 
embedding of the visible and hidden quarks into the same 
multiplet. The symmetry breaking pattern in this section 
is inspired by earlier models based on SU(3)w symme¬ 
tries, as described in [551 As described, we first 

consider all of the symmetries as global symmetries, as 
pNGB scenarios typically rely on gauging only some sub¬ 
group of the full global symmetry group this is in keeping 
with the standard structure of a pNGB scenario. We will 
return to the gauged subgroups later. We will also only 
consider the third generation fermions. 


TABLE III: U(l) Generators. 


Generator 

Diagonal Elements 

Al 

^(0,0,0,1,-1,0) 

A 2 

^(0,0,0,1,1,-2) 

A 3 

^(1,1,1,-1,-1,-1) 

A 4 

^(1,-1,0,0, 0,0) 

A 5 

^(1,1,-2, 0,0,0) 

As 

<ZB^(1, 1,1, 1,1,1) 

Ax 

qx ^{1,1,1,1,1,1) 


TABLE 11: Symmetries of the full matter content. 


Field 

SU(6) 

U(1)b 

SU(3)rv 

U(l)x 

Q 

6 

1 

3 

-1/3 

[/" X 2 

6 

-1 

1 

0 


6 

-1 

1 

1 

L 

1 

0 

3 

-1/3 


1 

0 

1 

1 

Hi 

1 

0 

3 

-1/3 

H 2 

1 

0 

3 

-1/3 

(j) X 3 

6 

1 

1 

0 


In Table [II| we show the field content and full symmetry 
group of the relevant matter in a UV completion of the 
SU( 3 )rv model detailed in Sec. |IV[ The fields Hi^2 and 
4 ‘i, 2,3 are scalars. H12 will eventually be responsible for 
breaking SU( 3 )w x U(l)x —t SU( 2 )w x U(l)v, however 
this aspect will not be relevant for proton decay. As we 
will see, on the other hand, (/'i,2,3 will be responsible for 
breaking U( 1 )bxSU( 6 ) —)■ U( 1 )b^, xSU( 3 ) where 
is baryon number for the visible sector quarks. 

The SM couplings arise from the following terms 

C D QhIu^ + QhIu^ + ^QHiH2D‘^ + jLHiH2E‘^ , 

(Al) 

and the neutrino masses from 

. (A 2 ) 

Let us first consider the breaking of the full color group. 
We will assume a scalar potential at the cutoff where 
^1,2,3 all obtain VEVs in the last three components such 
that they break SU(6) —?► SU( 3 ). We can determine the 
remaining U(l) symmetries by considering the full set of 
Abelian generators acting on the quarks. We may choose 
an arbitrary basis for the SU(6) generators and we thus 
choose to consider the basis of generators described in 

Table rrm 

When {(1)1,2,3) ^ 0 the generators Ai_2.3,b all appear to 
be broken. However only one diagonal combination of the 
final two are broken and the full set of broken generators 
is Al, A2, (A3 — As). The unbroken ones are A4, A5, (A3 -I- 
Xb) and Ax. The generator Ax is unbroken as (j) carries 
no charge under this symmetry. The generators A4, A5 
correspond to the full unbroken SU( 3 ) symmetry which 


TABLE IV: Symmetries of the fields below the scale of SU (6) 
breaking. All quarks are now written in lower case and the 
subscript T denotes a hidden sector quark. Representations 
under the spontaneously broken SU(3 )t are also shown. 


Field 

SU(3), 

SU(3)t 

U(1)sv 

SU(3)w 

U(l)z 

<7 

3 

1 

1 

3 

1/3 

X 2 

3 

1 

-1 

1 

-2/3 

d" 

3 

1 

-1 

1 

1/3 

<7t 

1 

3 

0 

3 

-1/3 

Ut X 2 

1 

3 

0 

1 

0 

<Et 

1 

3 

0 

1 

1 

L 

1 

1 

0 

3 

-1/3 


1 

1 

0 

1 

1 

Hi 

1 

1 

0 

3 

-1/3 

H 2 

1 

1 

0 

3 

-1/3 


we identify as the color symmetry of the SM. The hnal 
generator (A3 -I- Ab) is the most interesting. By studying 
the quark charges under the parent symmetry U(1 )b x 
SU(6) shown in Table [Tlj one finds that the remaining 
U(l)3+B = U(1)B.^, symmetry is simply baryon number 
symmetry for the SM quarks. 

Thus, to summarise, if in addition to the SU(6) sym¬ 
metry we assume a global U(1 )b symmetry acting on 
all of the SU(6) quark multiplets in the usual way, then 
when the scalars <('1,2,3 obtain vevs to break SU(6) down 
to the SM color group an additional U(l) symmetry will 
remain unbroken. This U(l) symmetry is simply the 
baryon number symmetry acting on the quarks of the 
SM and, as it is unbroken, this symmetry protects from 
proton decay as the proton is the lightest state carrying 
this baryon number. 

We may continue to study this theory to understand 
how the electroweak symmetry of the SM emerges. We 
will again choose to consider a rotated linear combina¬ 
tion of the U(l) symmetries U(1)2 =x-i-2/3Bv x U(1)bv- = 
U(l)x X U(1 )bv Up to overall choices in normalization 
this has not broken any of the global symmetries, thus 
the original U(1 )bv is still a symmetry of the theory, 
however this choice will be convenient for understand¬ 
ing the electroweak symmetry structure. At present the 
matter content and symmetry representations under the 


unbroken symmetries are given in Table IV 
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TABLE V: Symmetries of the fields below the scale of SU(6) 
breaking and after SU(3 )h/ x U(1)z —t SU(2)vc x U(l)y break¬ 
ing. Fields which may be projected out by higher dimension 
boundary conditions are shown in red. 


Field 

SU(3)c 

SU(3)t 

U(l)sv 

SU(2)m^ 

U(l)v 

4 

3 

1 

1 

2 

1/6 

u"" 

3 

1 

-1 

1 

-2/3 


3 

1 

-1 

1 

1/3 

N 

1 

3 

0 

1 

0 


1 

3 

0 

1 

0 

L 

1 

1 

0 

2 

-1/2 

E‘= 

1 

1 

0 

1 

1 

h\ 

1 

1 

0 

2 

-1/2 

h2 

1 

1 

0 

2 

-1/2 

Si 

1 

1 

0 

1 

0 

S2 

1 

1 

0 

1 

0 

Ls 

1 

1 

0 

1 

0 

Tl 

3 

1 

1 

1 

2/3 

Tr 

3 

1 

-1 

1 

-2/3 

qx 

1 

3 

0 

2 

-1/2 


1 

3 

0 

1 

0 


1 

3 

0 

1 

1 


The relevant interactions are now 

C, Z) qU^U2 1 2 

jqHiH2d^ + jqxHiH2d^ + jLHiH2E^{A3) 
and the neutrino masses will arise from couplings 

+ . (A4) 

The extended electroweak symmetry structure now mim¬ 
ics the set up found in more conventional composite Higgs 
models such as [sassHiz] and it is clear that a pNGB 
with uncolored top partners will arise in this scenario. 
No fields which obtain a vev from this point onwards 
are charged under the U(1)b global symmetry, thus it 
is clear that proton decay is avoided at the level of only 
global symmetries. However, for the sake of completeness 
we will now describe how the electroweak gauge group 
emerges. 

Both of Hi ^2 obtain vevs and thus break SU(3)w x 
U(l)z —SU(2)w X U(l)y, where the U(l)y symme¬ 
try is a diagonal combination of the SU(3)w generator 
proportional to (1,1,—2) and the E{l)z generator pro¬ 
portional to (1,1,1). The charges of the remaining fields 
under this final symmetry group is shown in Table |V] 
and the Hypercharge is identified as H = Z + Tg, where 
r8 = i/6(-i,-i,2) 

This demonstrates that the required global symmetry 
structure is compatible with the electroweak gauge group 
and a conserved global U(1 )b baryon-number symmetry 
in the IR. 


1. Gauging subgroups 

While it is clear that the desired global symmetry 
structure may emerge in the IR, in pNGB scenarios sub¬ 
groups of the global symmetries must be gauged to re¬ 
alize the SM gauge symmetry. With regard to the elec¬ 
troweak gauge symmetry only the SU(2)w C SU(3)vy 
will be gauged. If the full SU(3)w symmetry were gauged 
the Higgs boson would be eaten and could not emerge as 
a pNGB. The U(I)y subgroup must also be gauged to 
realize SM hypercharge. This symmetry arose as a lin¬ 
ear combination of a diagonal SU(6) generator, the orig¬ 
inal U(1)b symmetry, and the original U(l)x symmetry. 
Thus the realization of hypercharge is necessarily linked 
to the full SU(6) symmetry structure. Thus we will only 
gauge the SU(3) x SU(3) C SU(6) factors of the color 
sector, and set the gauge couplings equal by an exchange 
symmetry. This leaves all of the required U(l) genera¬ 
tors, including the A3 oc I/-\/6(1,1,1, —I, —1, —I) com¬ 
ponent of the full SU(6) symmetry, untouched. We may 
then gauge the linear combination of these symmetries 
and identify it with hypercharge. U(1)B^, also emerges 
as an unbroken global symmetry as a consequence of the 
full SU(6) symmetry. 

Considering the matter content that remains in the 
UV, shown in Table |Vj we see that with this matter con¬ 
tent hypercharge is anomalous. We will in fact start out 
with the ‘twisted multiplets’ of Sec.|I^ Taking these split 
multiplets removes precisely those fields in red in Ta¬ 
ble |V] and will leave hypercharge anomaly-free. Simply 
removing these components of the full symmetry multi¬ 
plets represents a hard breaking of the global symmetry, 
however enough of the symmetry remains to preserve the 
cancellation of quadratic divergences. A UV justification 
for the removal of these states may be found in extra- 
dimensional theories and projection by boundary condi¬ 
tions. Crucially, this symmetry breaking does not intro¬ 
duce any breaking of the U(1)B^, which prevents proton 
decay, thus it is innocuous from this perspective. 

To conclude, we see that a stable proton is consistent 
with the full global symmetry structure in the UV, and 
with the breaking required to achieve the desired IR spec¬ 
trum. Furthermore, the symmetry which stabilizes the 
proton is also consistent with the hard global symme¬ 
try breaking introduced whenever subgroups of the full 
global symmetry are gauged. 

2. Relating to GUTs 

Finally we wish to comment on the difference between 
the models discussed here and proton decay in Grand 
Unified Theories (GUTs). Since the hidden sector quarks 
mix with the SM neutrinos the hidden sector quarks effec¬ 
tively become leptons, thus the situation described here 
and the situation for GUTs is very analogous. The source 
of proton decay in GUT theories is neatly described in 
|94) . In GUTs if quarks and leptons live in the same 
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multiplets proton decay does not arise. The reason for 
this is that an effective baryon number is carried by the 
GUT gauge bosons, however this is still conserved and 
this alone cannot convert e.g. qqq I + qq. The proton 
decays through the combination of quarks and leptons 
living in the same multiplet and the fact that quarks and 
antiquarks live in the same multiplet, i.e. in the 10 of an 
SU(5) multiplet or a 16 of SO(IO). It is this collective 
symmetry breaking that allow proton decay. One way of 
seeing this is that exchange of a GUT gauge boson may 


mediate qq ^ql. 

Returning to the Natural Neutrinos scenario, this sec¬ 
ond condition is not met. Thus in principle if we were 
to gauge the full SU(6) symmetry then there would ex¬ 
ist heavy gauge bosons which (after mixing through the 
neutrino mass terms) do couple quarks and leptons. Pro¬ 
cesses such as ql ql would be mediated by such gauge 
bosons. However these gauge bosons would not also me¬ 
diate qq —> ql as quarks and anti-quarks live in separate 
mnltipets and a global U(1 )b symmetry is preserved. 
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